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X-Ray Observations of RX J0822-4300 and Puppis-A 
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Abstract. Based on observations with the X-ray observatories Chandra and XMM-Newton we present results from a detailed 
spectro-imaging and timing analysis of the central compact X-ray source RX J0822-4300 in the supernova remnant Puppis-A. 
The superior angular resolution of Chandra allows for the first time to pinpoint the point source nature of this object down to 
0.59 + 0.01 arcsec (FWHM) and to determine its position: RA=08''21"'57.40% Dec=-43°00' 16.69" (J2000) with sub-arcsecond 
accuracy. Spectral fits based on Chandra and XMM-Newton data provide a tight constraint on the emission properties of 
RX J0822-4300. Most of its X-ray emission seems to be of thermal origin. A model spectrum consisting of two blackbody 
components with Ti =^ 2.6 x 10* K, T2 - 5.0 x 10* K and ^1 ^ 3.3 km, ^2 - 0.75 km for the blackbody temperatures and 
the size of the projected emitting regions, respectively, provides the best model description of its spectrum. A search for X-ray 
pulsations from RX J0822-4300 revealed an interesting periodicity candidate which, if confirmed, does not support a scenario 
of steady spin-down. 
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1. Introduction 

For many years, it has been generally believed that all young 
neutron stars have similar properties as those observed in 
young rotation-powered pulsars, i.e. emitting strongly pulsed 
plus powerful plerionic radiation caused by non-thermal emis- 
sion processes in the neutron star's magnetosphere. Many re- 
cent observations of compact X-ray sources in supernova rem- 
nants (SNRs), however, suggest that this picture is incom- 
plete and no longer justified. Apart from appearing as rotation- 
powered pulsars, it has been shown that there are other man- 
ifestations of young neutron stars, e.g. with no radio coun- 
terpart identified. There are slowly rotating (P~ 6 - 12 s) 
compact objects which possibly have an ultra strong (B ~ 
10'"* - 10'^ G) magnetic field. These neutron stars are dubbed 
as magnetars, which include the anomalous X-ray pulsars 
(AXPs) and the soft gamma-ray repeaters (SGRs), depend- 
ing on whether bursts of strong y-ray emission is detected 
from them (e.g., Mereghetti 1998; Thompson 2000). The other 
class of objects are the "radio-quiet neutron stars" (e.g. Brazier 
& Johnston 1999). Most of them were identified by their 
high X-ray to optical flux ratios, others simply by their lo- 
cations near to the expansion centers of SNRs (e.g. Becker 
& Pavlov 2001, Kaspi et al. 2004), strongly suggesting that 
they are indeed the compact stellar remnants formed in the 
supernova events. The group of SNRs which are known to 
host a radio-quiet but X-ray bright central compact object is 
a slowly growing one. Thanks to more sensitive X-ray obser- 
vatories it currently includes Cas-A (Tananbaum 1999), the 
Vela-Jr. remnant (RX J0852.0-4622; Aschenbach 1998), RX 
J1713.7-3946 (Pfeffermann & Aschenbach 1996), RCW 103 



(Tuhoy & Garmire 1980), Puppis-A (Petre et al. 1982), PKS 
1209-51/52 (Helfand & Becker 1984) and Kes 79 (Seward et 
al. 2003; Gotthelf et al. 2005). 

The discovery of the X-ray point source, RX J0822-4300, 
in Puppis-A was initially made in one of the EINSTEIN 
HRI images of the SNR G260.4-3.4 (Petre et al. 1982). 
RX J0822-4300 appeared in this data as a fwcA X-ray feature. 
With ROSAT, it became strongly evident that RX J0822-4300 
is the compact stellar remnant which was formed in the SN 
event (Petre, Becker & Winkler 1996; hereafter PBW96), al- 
though the positional offset from the SNR's optical expansion 
center is 6.1 arcmin (cf. Winkler & Kirshner 1985, Winkler et 
al. 1988). The age of Puppis-A, estimated from the kinemat- 
ics of oxygen-rich filaments is ~ 3700 years. The remnant's 
kinematic distance estimated from an HI study of the inter- 
stellar medium along the line of sight towards Puppis-A is 
~ 2.2 + 0.3kpc (Reynoso et al. 1995; 2003). The space velocity 
of RX J0822-4300 required to travel to its observed position 
thus is ~ 1000 km/s. This is very high if compared with the 
mean proper motion velocity observed in ordinary field pul- 
sars but still comparable with what is observed in several of 
the other young supernova/pulsar associations (Manchester et 
al. 2005). 

RX J0822-4300 has not been detected as a radio pulsar. 
Limiting radio flux densities at 436 MHz, 660 MHz and 1520 
MHz are 1.5 mJy, 1.3 mJy and 0.3 mJy, respectively (Kaspi 
et al. 1996). For comparison, the limiting flux density in the 
Parkes Multi Beam Survey along the galactic plane was ^ 0.2 
mJy (Manchester et al. 2001) and the typical limiting sensitiv- 
ity in deep searches for young radio pulsars in SNRs is 50 /^Jy 
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(Camilo 2003). Gaensler, Bock & Stappers (2000) searched for 
a radio nebula around RX J0822-4300 with a resolution of 
27 X 17 arcsec. Their non-detection of any extended plerionic 
radio emission up to a scale of 30 arcmin prompted them to 
conclude that if RX J0822-4300 is a rotation-powered pulsar, 
e.g. with the radio beam not intersecting with the observer's 
line of sight, then it must be less powerful than other typical 
young radio pulsars located in SNRs. All young radio pulsars 
which are associated with a SNR have a spin-down power in 
excess of ~ 10^*^ ergs/s and are observed to power a X-ray/radio 
bright pulsar-wind nebula. 

RX J0822-4300 has no optical counterpart down to a lim- 
iting magnitude of B > 25.0 and R > 23.6 (PBW96). This limit 
yields an X-ray-to-optical flux ratio fxlfs > 5000 (PBW96). 
Together with the radio upper limits this rules out many types 
of X-ray sources as a likely counterpart of RX J0822-4300, 
except a neutron star. 

PB W96 fitted the ROS AT PSPC spectrum with a blackbody 
model and obtained a temperature of (3.2 + 1.2) x 10* K and 
a column density of (4.1 ± 0.2) x 10^' cm"^. The radius of the 
corresponding blackbody emitting area in their fits is only ~ 2 
km. Zavlin, Triimper, & Pavlov (1999; hereafter ZTP99) tested 
whether a hydrogen atmosphere model could bring this result in 
better agreement with the predictions of standard cooling mod- 
els. They fitted a temperature which is about half that found by 
PBW96 though with an increased radius of 10 km for the emit- 
ting area. However, as atmosphere models are seen not to be in 
agreement with the spectral fits from the cooling neutron stars 
Geminga, PSR B0656+14 and PSR B1055-52 (e.g. De Luca et 
al. 2005), the applicability of those models which in most cases 
use non-magnetic opacities only, is restricted. 

EINSTEIN and ROSAT data do not show any evidence for 
short or long term flux variations. Although a marginal detec- 
tion of X-ray pulses at a period of ~ 75.3 ms was claimed 
by Pavlov, Zavlin, & Trumper (1999; hereafter PZT99), it 
could not be confirmed so far (Pavlov et al. 2002; Becker & 
Aschenbach 2002). 

In order to put tighter constraints on the emission prop- 
erties of RX J0822-4300, various observations with the new 
generation X-ray satellites XMM-Newton and Chandra were 
targeted to it in the past few years. Making use of XMM- 
Newton's huge collecting power and high spectral resolution as 
well as of Chandra's sub-arcsecond angular resolution we have 
performed a sensitive broadband spectro-imaging analysis of 
RX J0822-4300 and its environment using all XMM-Newton 
and Chandra data taken from this source so far. This is the sub- 
ject of this paper which is organized as follows. In §2 we give 
a brief description of the relevant XMM-Newton and Chandra 
observations. In §3, we present the methods and results of our 
data analysis which are discussed in §4 in the context of a num- 
ber of physical models for the nature of RX J0822-4300. 

2. Observations 

In total, five observations have been targeted with XMM- 
Newton and Chandra on RX J0822-4300. All data have been 
taken between December 1999 and November 2001. We sum- 
marize the basic information of these observations in Table [l] 



and give a more detailed description in the following subsec- 
tions. 

2.1. XMM-Newton Observations 

Two of the five data sets reported here were obtained with 
the European Photon Imaging Camera (EPIC) aboard XMM- 
Newton (Jansen et al. 2001). EPIC consists of two Metal Oxide 
Semiconductor (MOSl/2) CCD detectors (Turner et al. 2001) 
of which half of the beam from two of the three X-ray tele- 
scopes is reflected to. The other two halves of the incoming 
photon beams are reflected to a grating spectrometer (RGS) 
(den Herder et al. 2001). The third of the three X-ray tele- 
scopes is dedicated to expose the EPIC-PN CCD detector 
solely (Struder et al. 2001). The April 2001 XMM-Newton ob- 
servation (hereafter XMMl) was taken with a total exposure 
time of ~ 28.8 ksec. The November 2001 observation (here- 
after XMM2) had an exposure time of ~ 24.3 ksec. The EPIC- 
PN CCD was operated in both observations in small-window 
mode with a thin filter to block optical stray light. This data pro- 
vide imaging, spectral and temporal information. All recorded 
events are time tagged with a temporal resolution of 5.7 ms. 
The MOS 1/2 CCDs were setup to operate in full- window mode 
with a medium filter in the April 2001 observation and a thick 
filter in the November 2001 observation. The MOS 1/2 cameras 
provide imaging, spectral and timing information, though the 
later with a temporal resolution of 2.6 s only. 

For both XMM-Newton observations, the satellite was 
pointed to RA=08*'21'"56^ and Dec=-43°00'19" [J2000]) 
which places RX J0822-4300 at the optical axis in the EPIC- 
PN CCDs. The raw data from the EPIC instruments were 
processed with version 6.0.0 of the XMM Science Analysis 
Software. Examining the raw data from the EPIC-PN CCD for 
both XMMl and XMM2, we did not find any timing anomaly 
observed in many of the XMM-Newton data sets (cf. Becker & 
Aschenbach 2002; Kirsch et al. 2004). This provides us with 
opportunities for an accurate timing analysis. We created fil- 
tered event files for the energy range 0.3 keV to 10 keV for all 
EPIC instruments. A small fraction of X-ray events might be 
split between CCD pixels. In order to correct for this effect only 
those events were accepted for which the corresponding X-ray 
generated pattern was between - 12 in MOS cameras and be- 
tween - 4 in the EPIC-PN camera ' . We further cleaned the 
data by accepting only the good times when sky background 
was low and removed all events potentially contaminated by 
bad pixels. The effective exposure times after data cleaning are 
summarized in column 7 of Tabled 

In order to correct for the non-uniformity across the detec- 
tor and the mirror vignetting, exposure map is needed to rescale 
all parts of the image to the same relative exposure. This is cre- 
ated by using XMMSAS task EEXPMAP 

RX J0822-4300 is located in a patchy SNR environment. 
This makes the extraction of its source and background spec- 
trum diflicult. In order to maximize the signal-to-noise ratio for 
RX J0822-4300, we extracted its source spectrum from circles 

' For a detailed description of EPIC event grade selection, please 
see the XMM-Newton Users' Handbook. 
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with 18 arcsec radii in both, the MOSl/2 and EPIC-PN cam- 
eras. About 70% of all point source events are located within 
the selection region. Annular regions with radii between 20-35 
arcsec, centered at RX J0822-4300, were used to extract the 
background spectra. The background corrected count rates are 
listed in column 8 of Table \l\ Response files were computed 
for all data sets by using the XMMSAS tasks RMFGEN and 
ARFGEN. 



2.2. Chandra Observations 

Three of the five data sets on RX J0822-4300 were taken 
with the Chandra satellite (e.g. Weisskopf 2004). One obser- 
vation was performed by using the Advanced CCD Imaging 
Spectrometer (ACIS; Burke et al. 1997) whereas the other two 
exposures were done by using the High Resolution Camera 
(HRC; Zombecket al. 1995; Murray et al. 1997). For the data 
reduction we used ClAO 3.0.2. 

The Chandra HRC data were taken on 1999 December 21 
and 2001 January 25 for HRC-1 and HRC-S, respectively, with 
RX J0822-4300 placed ~ 0.3 arcmin off'-axis. In order to de- 
termine the event positions accurately, we started the analysis 
with level- 1 event files and corrected for the tap-ringing distor- 
tion in the HRC event position reconstruction. Apart from this, 
we also performed the de-gap correction to the event files so 
as to compensate the systematic errors introduced in the event 
positions by the algorithm used to determine the centroid of the 
charge cloud exiting the rear micro-channel plate of the HRC. 
Furthermore, aspect offset was corrected for the event files. 

The ACIS observation was performed on 2000 January 1 
using the front-illuminated (FI) ACIS-S2 chip with a frame 
time of 0.84 s. RX J0822-4300 is located ~ 2.7 arcmin off^-axis 
in this ACIS observation. In order to correct for possible pileup 
effects, we started our analysis again with level- 1 files as those 
have preserved a number of source events which could have 
been misidentified as afterglows of cosmic ray events in the 
standard processing of level-2 data (cf. Davis 2002). The sub- 
arcsecond resolution of Chandra allows to extract the counts 
for the spectral analysis from a circle with radius 2.6 arcsec 
(encircled energy ~ 99% for on-axis point sources). This se- 
lection radius minimizes the contamination from the supernova 
background emission. An annular region with radii between ~ 
2.6-5.3 arcsec, centered at RX J0822-4300 was chosen to ex- 
tract the background spectrum. Response files were created us- 
ing the tools MKRMF and MKARF of CIAO. The background 
corrected HRC and ACIS count rates of RX J0822-4300 are 
given in column 8 of Tabled 

From the ACIS-S2 and XMM-Newton data we found that 
the energy of the central source peaks at ~ 1 .5 ke V. With the 
peak energy and for the ofF-axis angle of ~ 0.3 arcmin we ex- 
tracted the desired point spread function (PSF) model images 
from CALDB 2.26 standard library files (Fl) by interpolating 
within the energy and off-axis angle grids by using CIAO tool 
MKPSF. Exposure maps for the corresponding images were 
generated by the tool MKEXPMAP. 



3. Data Analysis 

3.1. Spatial Analysis 

Composite images of the supernova remnant Puppis-A and 
its central region around RX J0822-4300, as seen by the 
ROSAT HRI, by XMM-Newton's MOSl/2 CCDs and by the 
Chandra HRC-I, are shown in Figureffl RX J0822-4300 is lo- 
cated at the center of these images. From Figure^ it can be 
seen that the hardest X-ray emission in the remnant is mainly 
contributed by the central compact object RX J0822-4300. 
Apart from this, we observed two more hard X-ray sources 
in the XMM-Newton MOSl/2 images. Their locations as in- 
dicated by circles in Figurellt) are RA= 08''22'"26.70', Dec= 
-43°10'25.99" (J2000) for the source located in the south 
and RA= 08'>22'"24.25^ Dec= -42°58'00.82" (J2000) for the 
northern source. The location of the northern source is close 
to the region which was suggested by Winkler et al. (1989) to 
be a second supernova within Puppis-A. These authors have 
observed an unusual swirl-like structure in optical images and 
interpreted this as a possible second supernova remnant. The 
center coordinate of this structure is at about RA= 08''22'"39\ 
Dec= -42°59'41" (J2000). The left box in FigureUJ) illustrates 
the field of view in their observations. The angular separation 
between the northern hard X-ray source and the center of the 
swirl-like structure is ~ 3.2 arcmin. From the spectral analy- 
sis of optical filaments, Winkler et al. (1989) estimated that the 
kinematic age of the proposed second SNR is < 800 years. If 
the northern hard X-ray source is correlated with this structure 
and this age estimate is correct it would require a space veloc- 
ity > 2000 km/s (for an assumed distance of 2.2 kpc) in order 
to travel to its observed location. An association thus would 
be unlikely due to this high space velocity. For the southern 
hard X-ray source its correlation with Puppis-A is unspecified 
though most likely this is a background source. The photon 
statistics does not support a detailed spectral analysis for these 
two sources. 

The XMM-Newton MOSl/2 false color image (Figurellt)) 
demonstrates nicely that the south-western part of the remnant 
as well as the region near RX J0822-4300 comprise mainly 
hard X-ray photons. This is different from other parts of the 
image which consist of soft X-rays from the hot supernova 
ejecta. As Puppis-A is located at the edge of the Vela super- 
nova remnant (distance ~ 0.25 kpc), and is located behind it, 
we speculate that there is intervening absorbing material from 
Vela along the line of sight which absorbs most of the soft X- 
ray photons of the south-western part of Puppis-A. This view 
is supported by Figure|2lwhich shows a belt of absorbing mate- 
rial crossing the whole Puppis-A supernova remnant from the 
south-western to the north-eastern direction (Aschenbach 1994 
and discussion therein). 

The high resolution X-ray image from HRC-I allows for 
the first time to examine the spatial nature of RX J0822-4300 
with sub-arcsecond resolution. However, we have found that 
the full width half maximum (FWHM) of the point spread func- 
tion (PSF) (~ 0.4 arcsec) generated from the library files is 
narrower than expected. This can be ascribed to the fact that 
the PSF library files are derived by a ray-tracing program in- 
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stead of obtained directly from the calibration data. Due to this 
caveat, it is legitimate to fit the image with the convolution of 
a 2-dimensional Gaussian function instead of a delta function. 
The radial profile of RX J0822-4300 is depicted in Figure H 
The solid curve represents the best-fit Gaussian model with 
the modeled PSF at 1.5 keV as a convolution kernel. The best- 
fit results yield a FWHM of 0.59 + 0.01 arcsec which is very 
close to the expected width of the Chandra PSF. This result ap- 
pears as the first evidence for the point source emission nature 
of RX J0822-4300. Moreover, we were also able to narrow 
down the position of this compact object to the smallest region 
that has never been obtained before. The best-fit gives us a po- 
sition for RX J0822-4300 which is RA= 08*'21"57.40^ and 
Dec= -43°00' 16.69" (J2000). The statistical error of the posi- 
tion introduced by the image-fitting is found to be 0.01 arcsec 
(1 - cr). The predominant uncertainty of the source position 
is given by the finite width of the PSF (~ 0.5 arcsec) and the 
average pointing accuracy of the satellite (~ 0.6 arcsec). The 
same position is obtained from the analysis of the HRC-S data. 
The position and point source character of RX J0822-4300, de- 
duced by using Chandra, are in agreement with what we found 
in the XMM-Newton data (PSF ~ 5 arcsec FWHM). 

From observations of the neutral hydrogen surrounding 
RX J0822-4300, Reynoso et al. (2003) found a depression in 
the /121-cm line emission near to RX J0822-4300. According 
to their interpretation this structure could be connected to the 
compact stellar remnant because of its symmetric appearance 
as well as because of its alignment with the remnants opti- 
cal expansion center and the position of RX J0822-4300. In 
order to search whether there is an X-ray structure near to 
RX J0822-4300 which correlates with this radio structure we 
have overlaid the radio contours from Reynoso et al. (2003) 
on the XMM-Newton and Chandra HRC-I image (cf. Figure 
|5jl. No clear correlation between the radio and X-ray structures 
is seen, though the patchy supernova environment makes any 
conclusion uncertain. 

3.2. Spectral Analysis 

We estimated the eff'ects of pileup in both XMMl and XMM2 
data by using the XMMSAS task EPATPLOT. Our results 
showed that all the EPIC data were not affected by CCD pileup. 
Using the spectral parameters of RX J0822-4300 inferred 
from XMM-Newton, we estimated with the aid of PIMMS 
(version 3.6a) that the ACIS-S2 data are piled-up by a fraction 
of ~ 11%. We applied adequate correction by incorporating a 
pileup model in the spectral fitting (Davis 2001). Chandra data 
were also corrected for the degradation of quantum efficiency. 

In order to constrain the spectral parameters tightly, we fit- 
ted XMMl, XMM2, as well as the ACIS-S2 data simultane- 
ously. In order to obtain spectra from different observations and 
instruments with compatible significance, the energy channels 
were grouped dynamically with respect to the photon statistics 
in the analyzed data sets. For the MOSl/2 data of XMM2, we 
grouped the data to have at least 50 counts per bin. For the 
MOSl/2 data of XMMl as well as the ACIS-S2 data we ap- 
phed a grouping so as to have at least 100 counts per spectral 



bin. For the EPIC-PN data from XMMl and XMM2 a group- 
ing of 200 cts/bin was used. All spectral fits were performed in 
the 0.3 - 10 keV energy range by using XSPEC 11.3.1. 

Various model spectra like single blackbody, double black- 
body, power-law, combinations of blackbody and power-law, 
broken power-law, as well as thermal bremsstrahlung were fit- 
ted to the data. Independent of the fitted spectral models we 
found that the fits improve if the spectrum extracted from the 
Chandra data was not included in the analysis. Since the spec- 
tra extracted from XMM-Newton data supersedes the Chandra 
data in photon statistics, we excluded the later from the spectral 
analysis without loss of generality. The parameters of all fitted 
model spectra are summarized in Table 2. The quoted errors are 
conservative and are Icr for 2 parameters of interest for single 
component spectral models and for 3 parameters of interest for 
multi-component model. 

Fitting the spectral parameters of RX J0822-4300 as in- 
fen-ed from ROSAT PSPC data by PBW96 we found that these 
parameters (T = 3.2 x lO^K, Nh = 4.1 x lO^'cm"^^ = 2 km) 
yield no acceptable description of the XMM-Newton observed 
spectrum (xl - 21.17 for 467 dof). In general, spectral fitting 
with a single component blackbody or power-law model did 
not model the data beyond ~ 3 keV (cf. Table 2). Testing multi- 
component models we found that a two component blackbody 
with Nh = 4.54+JJ-4^ x lO-'cm^^, = 2.61+JJfg' x 10% 
T2 = 5.04+_°jl X 10% and Ri = 3.29!^;^^ ^j^^ 0.75!JJ;[2 
km for the blackbody temperatures and emitting areas, respec- 
tively, yields the best description of the observed spectrum. The 
reduced-;i'^ of this fit is 1.20 for 465 dof. We note that the ap- 
parent deviation of the reduced-;^'^ from one, indicating an ac- 
ceptable fit, is due to the fact that the data from different instru- 
ments and different epochs are modeled simultaneously. The 
benefit of combining all spectral data in simultaneous fits is the 
higher photon statistics and thus the ability to better discrimi- 
nate between competing model spectra. 

Figure |6l and Q shows the spectral fit for an absorbed dou- 
ble blackbody model and the corresponding spectral compo- 
nents respectively. In order to properly constraint the parame- 
ter space for the best-fitting model, we calculated the contour 
plots in the Ti - R\, T\ - Nh and T2 - Ri, T2 - Nh planes, 
respectively. These plots are depicted in Figure|S] For a consis- 
tence check we modeled the Chandra ACIS-S2 spectrum with 
the best fitting double blackbody model and found all parame- 
ters in agreement with those fitted for the XMM-data. 

Both Ri and R2 inferred from the double blackbody fit are 
inconsistent with the size of a canonical neutron star (i.e. R ~ 
10 km). It is therefore instructive to redo the fitting for this 
model with 7? 1 fixed at 10 km. This model still yields acceptable 
values of Nh = 6.38+[J-2i x lO^'cm^^, Ti = 1.87+° 02 ^ 10%, 
T2 = 4.58!°;°^ X 10%,' and R2 = 1.09!0;°4 km with;t-2 = 1.28 
(for 466 dof) only a slightly larger than leaving Ri uncon- 
strained. 

For a model combining a blackbody and a power-law, the 
goodness-of-fit is compatible with that for the double black- 
body model (xl = 1.21 for 465 dof). The inferred slope of 
the power-law component is F = 4.67^y yg. Although this is 
steeper than the photon-index, F = 1-3, observed for rotation- 
powered pulsars (cf. Becker & Triimper 1997), the model can- 
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not be rejected simply based on this as it is not a priori applica- 
ble for central compact objects. However, the column density 
is much higher than the expected level. When Nh is fixed to 
4 X 10^'cm"^, which is consistent with the values obtained by 
PBW96, ZTP99 and Winkler et al. (1981), it results in a param- 
eter set of r = 2.51!°;!^, T = 3.80!°;°^ x lO'-K, R = 1.59![;;[;^ 
km, though with a large of 1 .42 for 466 dof . 

It is necessary to examine whether a broken power law 
model can describe the spectra. This implies a purely non- 
thermal emission with spectral steepening at high energy af- 
ter an energy break which is due to the deficit of energetic 
emitting charged particles. From Table 2 it is obvious that the 
broken power law model does not yield any photon index that 
is consistent with that of a typical pulsar. We also fitted the 
data with a thermal bremsstrahlung model which physically 
implies that the central compact object would be surrounded 
by a hot plasma. From the normalization constant inferred from 
the spectral analysis, we can calculate the extent of the plasma. 
Following laria et al. (2001), we assume the bremsstrahlung 
normalization to be Ntremss - 3.02 x lO'^^N^V/4nD^, where 
D is the distance to the source in cm, A^^, is the electron density 
(cm^^), and V is the volume of the bremsstrahlung emitting 
region. Assuming A^^, is comparable with the average density 
~ Icm'^ of Puppis-A (Petre et al. 1982), the radius of the as- 
sumed spherical emitting region is estimated to be ~ 2.7 pc 
for an adopted distance of 2.2 kpc. This implies that the source 
should be extended (at a level of ~ 4.2 arcmin assuming a dis- 
tance of 2.2 kpc), in contradiction to the results from the spatial 
analysis. 

From both, the XMM-Newton MOSl/2 and the Chandra 
HRC-I images, some faint and diff'use hard X-ray emission 
around RX J0822-4300 seems to be present (cf. Figures ^ 
& |3}. Its nature can be determined by examining its spec- 
trum. We extracted the events in the X-ray filament near to 
RX J0822-4300 from the MOSl/2 cameras of XMMl from 
a 80 arcsecx30 arcsec box centered at RA=08^21™57.077', 
Dec=-43°01'15".42 (J2000). We found that its spectrum is 
consistent with an absorbed non-equlibrium ionization colli- 
sional plasma model (XSPEC model: VNEI) with goodness- 
of-fit of - 151.08 for 142 dof. The energy spectrum as 
fitted to this model spectrum is displayed in Figure |5] Line 
emission is easily recognized in this plot. The most obvious 
feature is the O VII and O VIII line complex at 0.662 keV 
and 0.65 1 keV respectively. Parameters inferred from the best- 
fitted model are the column density Nh = 3.70^g [2 x 10^'cm"^, 
the plasma temperature T = 1.62^^ }^j x 10''K, the ionization 
timescale r = 2.33^g }j x 10"s cm"-', and the metal abundances 
with respect to the solar values (O: 17.83+J-^^, Ne: 4.00+[]-^^, 
Si: 2.25+j-^^, S: 6.58+^-^^, Fe: 2.24+°-2'^) (quoted errors are Icr 
for 2 parameters of interest). These parameters imply a rela- 
tive abundance ratio 0:Fe to be about 6-9 times its solar 
value which strongly suggest an enhancement of oxygen in 
Puppis-A. This is in agreement with the conclusion drawn by 
Canizares & Winkler (1981). For the other elements, includ- 
ing H, He, C, Mg, Ar, Ca and Ni, we do not find any sign of 
enhancement and their abundances are in agreement with the 
solar values. We have performed the spectral fitting with dif- 
ferent selected backgrounds. Provide that the backgrounds are 



selected from low count regions, all the best-fitted values are 
found to be within the quote Icr errors above. Since the rem- 
nant environment is patchy and inhomogeneous, abundance ra- 
tios from different regions are not expected to be comparable. 
A detailed modeling of the variation of chemical abundance is 
beyond the scope in this paper A further detailed analysis of 
Puppis-A is in preparation and will be published elsewhere. 

Although the rim emission appears to be a part of the 
structure of Puppis-A, we also exam whether there is any 
non-thermal contribution in the emission by adding a Crab- 
like power-law component (i.e. with photon index of 2) in the 
spectral fit. The additional component does not improve the 
goodness-of-fit at all (x^ - 150.95 for 141 dof). A 3cr upper 
bound of the power-law model normalization is estimated to be 
6 X 10"^ photons keV"' cm"^ s"'. This implies a limiting flux 
of non-thermal plerionic emission, if any, to be 3.06 x 10"'^ 
ergs cm"^ s"' and 2.88 x lO"'-' ergs cm"^ s"' in 0.1-2.4 keV 
and 0.5-10 keV respectively. 

3.3. Timing Analysis 

3.3.1 . Search for long-term variabilities 

To check whether the energy fluxes measured from 
RX J0822-4300 by XMM-Newton and Chandra ai-e consis- 
tent with each other or whether there are significant long-term 
deviations observed in the different data sets we computed 
the flux for the best fitting double blackbody model from all 
available data. In order to compare the XMM-Newton and 
Chandra results with the existing ROSAT flux we restricted 
this computation to the energy range 0.1 - 2.4 keV. As 
shown in Table |3] all observed energy fluxes, from ROSAT 
to XMM-Newton are consistent with a constant energy flux 
of /,(0.1-2.4keV) ~ 3 x lO^'^ ergs cm^^ s"'. The same 
conclusion can be drawn from comparing the broadband fluxes 
from Chandra, XMMl and XMM2. The observed flux of these 
three observations are found to be ~ 4 x 10 ergs cm"^ s ' in 
the range of 0.5-10 keV. 

3.3.2. Search for coherent pulsations 

Although the lack of long-term variability and the spectral 
analysis eliminates some highly improbable models and hence 
helps us to put constraints on the properties of the central com- 
pact object, the most strong argument that this object is indeed 
a neutron star would come from the detection of X-ray pulsa- 
tions. Since the small-window mode was setup for the EPIC- 
PN camera in both XMMl and XMM2 observations, the 5.7 
ms temporal resolution of this data is sufficient to search for 
coherent short-term pulsations. 

The arrival times in both event files were barycentric cor- 
rected using the XMMSAS task BARYCEN. In order to min- 
imize the systematic errors induced in the barycentric correc- 
tion, we use the position inferred from the Chandra HRC-I im- 
age fitting (cf. §3.1) for correcting the arrival times in both 
data sets. The initial period searches were performed by ap- 
plying a fast Fourier transformation (FFTs) on both sets of 
photon arrival times separately. The advantage of having mul- 
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tiple data sets from different epochs supports to cross-check 
any periodicity candidates easily and prevents wrong identi- 
fications. FFTs were calculated for each time series with 20 
different binnings. Searches in the frequency domain were lim- 
ited at 0.01 Hz < / < 100 Hz. Promising frequency peaks 
appearing in the power spectra of both observations were cross- 
correlated and selected for subsequent searches using standard 
epoch-folding analysis. 

PZT99 reported the detection of a periodic signal with 
P = 75.2797300 ms and P = 1 .49 x lO^'^s s"' in their ROSAT 
analysis. We searched for coherent pulsations in a period range 
extrapolated to these spin parameters in XMM-Newton data. 
The photon statistics of this data, compared to the ROSAT data, 
is about a factor 25 higher. A similar analysis as reported in the 
present work was already performed by Becker & Aschenbach 
(2002) who could not confirm the existence of a pulsed signal 
in an extrapolated period range, neither in the ROSAT nor in 
the XMM-Newton data. We herewith fully confirm their result 
in our independent analysis. 

Taking peaks in the power spectra as initial candidates, we 
made a more detailed search using Z^, test where m is the num- 
bers of harmonics (Buccheri et al. 1983). We have detected 
periodicities of P ^ 0.218749 + 0.000001 s in XMMl (MJD 
52014.4634390827268 days) ^ and P = 0.222527 + 0.000002 
s in XMM2 (MJD 52221.8938398198225 days) which both 
have very similar properties. The quoted uncertainties indicate 
the Fourier resolution P^ jT in the corresponding observation, 
where T is the time span in the data set. Using the //-test (De 
Jager, Swanepoel, & Raubenheimer 1989), we found that H 
is maximized for the first harmonic. The calculated for the 
detected signals in XMMl and XMM2 are 28.10 and 28.92, 
respectively. The nominal probabilities for the identification of 
these signals by chance are 8 x 10"^ and 5 x 10"^, respectively. 
The pulse profiles and Zj distribution is given in Figure [TO] 
Both lightcurves are similar to each other and share the same si- 
nusoidal characteristics. Following Becker & Triimper (1999), 
we calculated the pulsed fraction of this signal by the bootstrap 
method proposed by Swanepoel, de Beer, & Loots (1996) and 
obtained P/ = 5 + 1 % in both, the XMMl and XMM2 data sets. 
The period time derivative calculated from the separation of the 
epochs of the two data sets is P = (2.112 + 0.002) x 10"'°s s"^ 

Taking the number of 10^ trials into account, the proba- 
bilities for finding these signals by chance is (5 - 8) x 10"^. 
However, given the similarity in period, pulse shape, signal 
strength and pulse fraction together with the detection of the 
signal in two independent XMM-Newton data sets makes this 
signal a very promising candidate to test and search for in fu- 
ture observations. 

In order to minimize the probability of a false detection 
we have investigated the possibility that the signals are in- 
duced from the readout processes in the detector CCDs or other 
cyclic processes operating during data acquisition. For this we 
have extracted events from Puppis-A in both, the XMMl and 
XMM2 data sets from a location near to RX J0822-4300, 
i.e. from CCD columns located at the same level in the readout 



^ The mean epoch of the observation in TDB at the solar system 
barycenter (SSB) 



direction as RX J0822-4300. The same procedure of timing 
analysis was applied on these events as applied to the events 
from RX J0822-4300. However, we did not detect any cyclic 
signal at a period near to ~ 0.22 s. We therefore can rule out 
that the detected pulsations are due to periodic systematics in 
the on-board data processing during data acquisition. 

To further cross-check this periodicity detection, we uti- 
lized an independent data set from Chandra HRC-S. Since the 
observation with HRC-S was performed in "imaging" mode 
where the outer segments of the micro-channel plate were dis- 
abled, the total count rate is below the telemetry saturation 
limit, so that all events can be assigned with accurate time 
and the HRC timing anomalies are minimized. This enables 
us to perform an accurate timing analysis on this data set, 
though the photon statistics is a factor of ~ 6 lower than in 
the XMM-Newton data. The event file of HRC-S was firstly 
barycentric corrected (with the position given in the HRC-I 
image fit) by CIAO tool AXBARY. With the P and P esti- 
mated from XMM-Newton data, we extrapolated the period to 
the epoch of HRC-S observation as an initial starting point. 
A detailed search around this period gives a promising can- 
didate at P = 0.217303 + 0.000002 s in the Chandra HRC-S 
data (MJD 51934.6266560833901 days). The //-test indicates 
the highest probability for Z^ = 45.94 wich yields a nominal 
chance probability of 3 x 10"^. This is not strong enough to con- 
clude a significant signal in the HRC-S data though we point 
out that there are only ~ 6000 counts available for this test. For 
a 5% pulsed fraction as indicated in the XMM-Newton data 
only ~ 300 counts would contribute to the pulsed component. 
The low significance of the signal found in the Chandra data 
thus would be in line with a low significance of the periodic 
signal. Archival ROSAT and ASCA data of RX J0822-4300 
are of small photon statistics so that we did not include this 
data in order to search for a pulsed signal near to 0.22 s. 

4. Discussion & Conclusion 

The lack of any detectable long term variability together with 
the high X-ray/optical flux ratio and the observed spectral char- 
acteristics makes it very unlikely that RX J0822-4300 is some- 
thing else but the compact stellar remnant formed in the core 
collapsed supernova which was left behind Puppis-A. Our 
image analysis shows that the compact object is the hardest 
X-ray source in the 30 arcmin central region of Puppis-A. 
Chandra and XMM-Newton data do not show any extended 
X-ray emission which could be plerionic emission powered by 
the compact remnant, though this is difficult to quantify given 
the patchy environment in which RX J0822-4300 is located. 

From the spectral fitting, we found that the point-source 
spectrum is compatible with a two component blackbody 
model. The best-fit model yields A^^ = 4.54+[}^9 x 102'cm~^ 
temperatures of Ti = 2.61+[5-^°xlO*K and = 5.04+°-2fjXlO*K 
for the projected blackbody emitting areas with radii Pi - 
3.29;^q '4 km and P2 = 0.75;^g [5 km, respectively. Compared 
to standard cooling curves (e.g. Yakovlev et al. 2004) Ti is a 
little higher than ~ 1.8 x 10'' K which would be expected for 
a 1.35 Mq neutron star with a stiff equation of state (Prakash 
et al. 1988). Pi is not quite consistent with the size of a typ- 
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ical neutron star. This was already found by ZTP99 who at- 
tempted to obtain a set of reasonable neutron star parameters 
by modeling the ROSAT data with spectral models which take 
the presence of a hydrogen atmosphere into account, though 
the applicability of these models are restricted. 

In the present work we have shown that a parameter set 
which is consistent with a standard cooUng neutron star model 
can be obtained by fixing Ri at 10 km. This yields a col- 
umn density of 6.38;|;q j3 x lO^^cm"^, temperatures of Ji = 
1.87;^^;^^ X lO^K and 7^2 = 'i-^Stoij ^ l^^^ for a projected 
blackbody emitting area with radius R2 = 1 .09+[;{;^ km. 

The results of our analysis suggest that the low tempera- 
ture component is emitted from a large fraction of the neutron 
star surface, while the high temperature component is emitted 
from a much smaller and hotter region. This double blackbody 
model, though, could be a two-step adaption for a wider tem- 
perature distribution which centrally peaks. Such a hot spot on 
the neutron star surface can be produced by several mecha- 
nisms. One of them is the bombardment of the polar cap re- 
gions by energetic particles accelerated in the magnetosphere 
backwards to the neutron star surface (Cheng, Ho & Ruderman 
1986; Cheng & Zhang 1999). Another way to produce a hot 
spot on the stellar surface is by anisotropic heat transport (Page 
1995). Since the heat conduction inside a neutron star is much 
more efficient in the direction along the magnetic field lines 
than that in the perpendicular direction, a complete model of 
coohng magnetic neutron star should lead to an anisotropic 
heat flow and hence produce hot spots on the stellar surface. 
In this scenario, one should expect the emission to be pulsed 
at the rotation period of the star as the hot spot goes across the 
line of sight. A pulsed X-ray flux as revealed by the putative pe- 
riodic signals seen in XMMl and XMM2 thus would support 
this scenario. 

The pulsed fraction of the putative periodic signal is 5 ± 1% 
in XMMl and XMM2. The emission from young pulsars like 
the Crab is compatible with being 100% pulsed (Tennant et 
al. 2001), whereas the fraction of pulsed photons is ~ 7% 
for the Vela pulsar and ~ 20 - 40% in many of the other X- 
ray detected pulsars (e.g. Becker & Pavlov 2001; Becker & 
Aschenbach 2002; Kaspi et al. 2004). A low pulsed fraction, 
however, is not unexpected though. When the general rela- 
tivistic effect is taken into account (Page 1995; Hui & Cheng 
2004), the pulsations are found to be strongly suppressed and 
the pulsed fraction is highly dependent on the mass to radius 
ratio of the star, the orientation of the hot spot and the view- 
ing angle geometry. This is due to the fact that the gravitational 
bending of light will make more than half of the stellar surface 
become visible at any instant and hence the contribution of the 
hot spot will be hampered. If the orientation of the hot spot is 
deviated from that of an orthogonal rotator and/or the star has a 
high mass to radius ratio, then a very low amplitude pulsations 
is expected, which makes the periodicity search difficult. 

The P deduced for the candidate periodicity would be 
among the largest spin-down rates in the neutron star popu- 
lation. The largest known P was inferred from SGR 1806-20, 
/> = (8 - 47) X lO-'^s s-', (Kouveliotou et al. 1998; Woods et 
al. 2002). If the identifications of P and P are correct, it implies 
a non-steady spin-down behavior of RX J0822-4300. This 



phenomenon is not unobserved. There are two SGRs (SGR 
1806-20 and SGR 1900+14) which show large changes in the 
spin-down torque up to a factor of ~ 4 (Woods et al. 2002). 
Moreover, deviations from a steady spin-down were also ob- 
served in the other radio-quiet neutron stars such as in SNR 
PKS 1209-51/52 (Zavhn, Pavlov, & Sanwal 2004). However, 
the rotational dynamics cannot be determined without ambigu- 
ity here and further observations are needed to confirm and/or 
refine this putative periodicity. 

An alternative proposal to explain the origin of X-rays 
from radio quiet compact objects in supernova remnants is ac- 
cretion onto a neutron star (e.g. Pavlov et al. 2000). In this 
scenario, the observed luminosity L is powered by an accre- 
tion rate of M = L/((c^). ( is the accretion efficiency which 
is expressed as ^ = O.lMi ^R'^^. Equating the expression of 
M wiffi Bondi formula (i.e. M = AnG^M^pv'^), we can ex- 
press the relation of the circumstellar baryon density as n = 
8 X 10^VjQQ(0.2/^)Mj^4L33 cm"^, where vioo is ffie velocity of 
ffie neutron star in ffie unit of 100 km s~' and L33 is ffie luminos- 
ity in the unit of 10^^ erg s ^ The offset of ~ 6.1 arcmin from 
the optical expansion center, estimated distance (~ 2.2 kpc) and 
estimated age (~ 3700 years) suggest a transverse velocity of 
~ 985 km s"^ Even we take L33 = 1, which is lower than ffie 
value inferred from the acceptable spectral fit L33 ~ 5, ffie ex- 
pression of n impUes a density of about 6 orders of magnitude 
higher ffian ffie expected value found by Petre et al. (1982). 

It is therefore safe to reject the scenario that the observed 
X-rays are powered by accretion from circumstellar matter. The 
stringent optical limit also rules out ffie possibility ffiat the 
accretion is from a massive companion. However, we cannot 
completely exclude ffie possibilities that the central object is 
accreting from a very close dwarf star or from a fossil disk 
(van Paradijs et al. 1995) which remained after the supernova 
explosion. For the first possibility, even though such a compact 
system is unlikely to remain bound in the disruption of ffie high 
mass progenitor (> 25 M© Canizares & Winkler 198 1), a deeper 
search for optical/lR counterpart can still be beneficial. 

As has been shown in this work, ffie properties of the com- 
pact object in Puppis-A are similar to those of other radio 
quiet neutron stars in many aspects (e.g. Cas A; Pavlov et 
al. 2000; Chakrabarty et al. 2001; Murray et al. 2002). Even 
though the nature of this class of object has not yet been com- 
pletely resolved, ffieir existence has already revolutionized ffie 
conventional notion of neutron stars and ffieir environment. 
Since it is easier to detect and identify active radio/X-ray pul- 
sars ffian the radio-quiet neutron stars which are only observ- 
able in X-ray and located in a patchy X-ray bright supernova 
remnant, it is plausible that they are more common than canon- 
ical pulsars. It is not unlikely to assume ffiat ffie relatively small 
number of detected radio-quiet neutron stars is due to observa- 
tional selection efifects. In this sense, identifying the nature of 
this class, including RX J0822-4300, is very important. 
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Table 1. Details of the XMM-Newton and Chandra observa- 
tions of RXJ0822-4300 



Telescope 



Instrument Instrument Mode 



Filter Obs.ID Start Date Effective Net Rate 

Exposure (cts s"') 



Chandra HRC-I 

Chandra ACIS-S 2 

Chandra HRC-S 

XMM-Newton MOSl 

XMM-Newton MOS2 

XMM-Newton PN 

XMM-Newton MOSl 

XMM-Newton M0S2 

XMM-Newton PN 



Default UVIS 749 1999-12-21 16 ks 0.21±0.01 

Timed Exposure OBF 750 2000-01-01 11 ks 0.71±0.01 

Imaging UVIS 1851 2001-01-25 19 ks 0.32±0.01 

PrimeFullWindow Medium 0113020101 2001-04-15 7.3 ks 0.50±0.01 

PrimeFulIWindow Medium 0113020101 2001-04-15 7.5 ks 0.52±0.01 

PrimeSmallWindow Thin 0113020101 2001-04-15 16 ks 1.62±0.01 

PrimeFullWindow Thick 0113020301 2001-11-08 7.5 ks 0.45±0.01 

PrimeFullWindow Thick 0113020301 2001-11-08 6.2 ks 0.46±0.01 

PrimeSmallWindow Thin 0113020301 2001-11-08 16 ks 1.63±0.01 



Table 2. Spectral Parameters inferred from data from XMMl and XMM2. 



Parameter 


BB 


BB+BB 


BB+BB" 


PL 


BB+PL 


BB+PL* 


BB+PL^ 


BKPL 


BREMSS 


Af/;(10^'cm"^) 


2 67+°"' 

-0.09 


4 54+0.49 
•- ^ -0.43 


6 38+°-2i 

-0.13 


9 26+° i'i 

^' -0.17 


8 57+°-** 

-0.35 


4.00 


8.02 
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-0.03 


2 61+°-'° 


1 87+°°2 




4 41+0.05 

• -^-0.07 


3 80+° °* 


1.87 




57+0.18 

^•-"-0.06 


(lO'^Kl 




5 04+°-2^ 


4 58+°"? 

^•-^"-0.07 
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89+°"* 

-0.08 


1 cq+o.06 

i-'^_0.06 
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(\ 7«+0.12 
















"bremss cm ; 


















1 98+002 


Py^oA^iAtev (ergscm-2 s"') 


5.82 X 10-12 


9.67 X 10-12 


2.04 X 10-11 


1.77 X 10-' 


2.60 X 10-' 


1.18 X 10-11 


2.18 X 10-1° 


3.45 X 10-1° 


2.40 X 10-11 


^Xo.5_iotev (ergs cm-2 s"') 


6.05 X 10-12 


8.94 X 10-12 


1.48 X 10-11 


4.46 X 10-11 


3.88 X 10-11 


8.28 X 10-12 


2.69 X 10-11 


2.48 X 10-11 


1.27 X 10-11 


Reduced 


1.82 


1.20 


1.28 


1.91 


1.21 


1.42 


3.77 


1.22 


1.20 


D.O.F 


467 


465 


466 


467 


465 


466 


467 


465 


467 



The emitting areas inferred from the blackbody fits are calculated for the assumed distance of 2.2 kpc. 

Fxo i_24fav and Fxo^-mkev are unabsorbed fluxes in the energy range 0.1 - 2.4 keV and 0.5 - 10 keV, respectively. 

^ R\ fixed at 10 km radius. 

i' A^Hfixedat4x lO^'cm-^. 

" Ri and Ti fixed at 10 km and 1.87 x 10^ K respectively. 
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Table 3. Observed fluxes inferred from observations at different epochs. 



11 



Observations Fx (0.1-2.4 keV) Fx (0.5-10 keV) Start Date 
ergs cm"^ s"' ergs cm"^ s"' 



XMMl 


3-37!o:26 X 10- 


-12 


4-19!a36 X 10- 


-12 


2001-04-15 


XMM2 


3.3S-!_°-l° X 10- 


-12 


4.17!°;«x 10- 


-12 


2001-11-08 


Chandra ACIS 


3.34;0;^o X 10- 


-12 


4.18!^;^* X 10- 


-12 


2000-01-01 


ROSAT PSPC " 


~ 3 X 10^'^ 








1991-04-16 



" Flux observed by ROSAT according to PBW96 




Fig. 1. a. Composite ROSAT HRI image of the Puppis-A supernova remnant. The blue ring indicates the 30 arcmin central 
region which has been observed by XMM-Newton in April and November 2001. b. XMM-Newton MOSl/2 false color image 
of the inner 30 arcmin central region of Puppis-A (red: 0.3 - 0.75 keV, green: 0.75 - 2 keV and blue: 2-10 keV). The central 
source is RX J0822-4300. The inset shows the squared region as observed by the Chandra HRC-I. It is interesting to note that 
the region around RX J0822-4300 comprise mainly hard X-ray photons. The left white box indicates the region of the swirl-like 
structure interpreted by Winkler et al. (1989) as a second supernova in Puppis-A. The location of the northern and southern hard 
X-ray point sources are indicated by circles. The binning factors in the XMM and Chandra images are 4 arcsec and 0.6 arcsec, 
respectively. Adaptive smoothing with a Gaussian kernel of cr < 1 pixel has been applied to the XMM-Newton and Chandra 
images. Top is north and left is east. 
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Fig. 2. False color image of Puppis-A and parts of the Vela supernova remnant as seen in the ROS AT all-sky survey. Puppis-A 
appears solely in blue color (harder X-rays) while the soft emission from the Vela supernova remnant is represented by the red and 
yellow colors. The image clearly demonstrates the existence of the absorption belt which crosses Puppis-A from the south-west 
to north-east direction and which is associated with rim emission from the Vela supernova remnant. 




Fig. 3. False color image (red: 0.3 - 0.75 keV, green: 0.75 - 2 keV and blue: 2-10 keV) of the 4.4' x 4.4' region around RX 
J0822-4300 as seen by XMM-Newton's EPIC-PN detector during the observations in April 2001 (a.) and November 2001 (b.). 
The rim emission from Puppis-A near to the location of RX J0822-4300 is visible very well in a. Top is north and left is east. 
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Fig. 4. Chandra HRC-I brightness profile for RX J0822-4300. The solid curve represents the best-fit Gaussian model with the 
modeled PSF at 1 .5 keV as a convolution kernel. 
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Fig. 5. The contours depict a double-lobed feature of reduced /121-cm line emission (from Reynoso et al. 2003) superimposed on 
fluxed image of XMM-Newton MOSl/2 (upper) and Chandra HRC-I (lower) in gray scale. The cross shows the optical expansion 
center of Puppis-A as calculated by Winkler et al. (1988). 
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Fig. 6. Energy spectrum of RX J0822 -4300 as observed with the EPIC-PN (upper spectra) and MOS 1 /2 detectors (lower spectra) 
and simultaneously fitted to an absorbed two component blackbody model (upper panel) and contribution to the^^ fit statistic 
{lower panel). 




Fig. 7. Components and combined model of double blackbody fitted to the spectral data of RX J0822-4300. 
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Fig. 8. Icr, 2o" and 3cr confidence contours for the double blackbody fit to the X-ray spectrum of RX J0822-4300. 
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Fig. 9. Energy spectrum of the rim emission underneath RX J0822-4300 as observed in April 2001 with the EPIC-MOS 1/2 
detector and simultaneously fitted to an absorbed non-equiUbrium ionization coUisional plasma model (upper panel) and contri- 
bution to the;^^ fit statistic {lower panel). 
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